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a practical  technique  to  form  useful  wrought  products  from  materials  difficult  or 
impossible  to  work  by  other  processes,  and  it  offers  the  means  to  develope  superior 
useful  properties.  This  program  of  research  is  concerned  wich  the  response  of  a 
variety  of  materials  to  hydrostatic  extrusion.  A major  objective  j.s  to  relate  the 
microstructure  and  mechanical  properties  of  extruded  materials  to  the  important  pro- 
cessing variables,  extrusion  ratio,  extrusion  pressure,  temperature,  and  extrusion 
rate.  The  experiments  fall  into  two  categories;  1.  A detailed,  systematic  investi- 
gation of  the  influence  of  the  processing  variables  on  the  properties  of  several 
fairly  simple  materials  (pure  Fe,  Ni,  Mg,  and  Ti:  and  Cu-30Zn  or  Cu-lO^n)  2.  The 
utilization  of  hydrostatic  extrusion  to  optimize  the  pro] erties  of  more  complex 
materials  chosen  for  their  technological  utility  (Mg-Li-B  for  high  specific  stifiness; 
Fe-C  for  high  strength,  ductility;  T-D  Nichrome  for  high  temperature  strength; 

Al-Fe,  a high  conductivity,  low  density  material;  NbTi  a super  conducting  alloy). 

The  theoretical  portion  of  the  program  is  devoted  to  an  understanding  of  extrusion. 
Elastic-plastic  finite-element  programs  are  being  considered  as  a means  oi  analyzing 
the  extrusion  process.  Such  studies  are  directed  towards  a complete  stress  analysis 
of  extrusion  and  should  lead  to  a better  elucidation  of  the  factors  influencing  the 
onset  of  undesirable  defects  during  extrusion. 
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INFLUENCE  OF  HIGH  HYDROSTATIC  PRESSURE  EXTRUSION 

ON 

MECHANICAL  BEHAVIOR  OF  MATERIALS 


Introduction  and  Summary 

Hydrostatic  extrusion  is  a promising  new  metal  working  technique.  In 
hydrostatic  extrusion  a large  hydrostatic  pressure  is  superimposed  on  the 
stress  field  that  normally  develops  during  extrusion.  As  a result,  materials 
difficult  or  impossible  to  extrude  by  the  conventional  extrusion  process  can 
be  extruded  hydrostatically.  Further,  there  is  evidence  that  materials  formed 
by  hydrostatic  extrusion  exhibit  superior  mechanical  properties  compared  to 
similar  materials  worked  by  conventional  techniques.  At  the  present  time 
hydrostatic  extrusion  is  not  being  applied  extensively  to  technological 
applications  in  this  country  (considerably  more  is  being  done  in  Europe 
and  Japan).  A pilot  plant  to  extrude  hydrostatically  copper  into  wire  has 
been  started  by  Western  Electric  and  this  particular  operation  appears  very 
promising.  Surprisingly  the  basic  understanding  of  the  hydrostatic  extrusion 
process,  both  from  the  theoretical  mechanics  point  of  view  and  from  the 
standpoint  of  materials  science,  remains  in  quite  a prirxtrve  state. 

The  present  program  is  centered  on  an  experimental  and  analytical  study 
of  hydrostatic  extrusion.  The  hydrostatic  extrusion  press  and  auxiliary 
equipment  arrived  in  July  1975  and  was  made  operational  by  late  August  1975. 
The  experimental  aspects  are  described  in  the  first  section  of  this  progress 
report  (pages  2-11).  In  the  second  section  (pages  12-31)  we  describe  the 
basis  for  an  elastic-plastic  code  as  a means  of  assessing  metal  forming 
operations.  This  analytical  approach, is  unique  and  promises  to  be  a very 
powerful  tool  in  analyzing  the  extrusion  process. 
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I . Experimental  Aspects  of  Hydrostatic  Extrusion  Program 

(Robert  Whalen,  Research  Associate  and  Larry  Eiselstein,  Research  Assistant) 


After  a series  of  delays  caused  by  a variety  of  factors  such  as  machining  and 
supply  problems,  modifications  to  the  design,  and  an  overly  optimistic  original 
promise  date,  the  hydrostatic  extrusion  machine  was  finally  delivered  to  Stanford 
University  in  June  of  1975.  The  major  thrust  of  our  activities  at  Stanford  since 
that  date  has  been  on  the  installation  and  investigatory  testing  of  the  extruder. 

The  machine’s  special  power  requirements  necessitated  the  re-wiring  of  a portion 
of  the  John  Blume  Structures  Laboratory  in  which  the  extruder  is  located.  An 
investigation  of  alternatives  in  this  area  yielded  two  options.  The  first  of  these 
was  to  run  an  80  Amp  capacity,  three  phase  line  to  the  hydrostatic  extruder  and 
operate  as  a 220  volt  system.  The  second  option  was  to  utilize  a step-up  transformer 


which  would  provide  440  volts  at  the  main  box  and  require  the  installation  of  a 40  Amp, 
three  line.  Fortunately,  we  discovered  a surplus  power  transformer  on  campus  and 


this  brought  estimates  for  each  of  the  two  possibilities  to  approximately  the  same 
cost.  There  was  another  consideration  which  influenced  our  thinking  and  finally 
proved  persuasive.  *he  extruder  already  shares  the  Structures  Lab  with  other  heavy 
power  consuming  machinery  and  computers.  As  the  laboratory  expands  there  should  be 
a definite  advantage  to  maintaining  a distinct  power  system  which  would  be  independent  | 
of  power  demands  and  fluctuations  caused  by  these  competing  users.  Therefore  the 
decision  was  made  to  install  the  surplus  transformer  at  the  main  line  and  run  the 
extrusion  project  on  440  volts. 

Location  of  the  machine  in  the  Structures  Lab  is  adequate  for  our  purposes. 

The  site  provides  sufficient  work  space  for  our  purposes  and  the  Structures  Lab 
is  close  to  other  essential  areas  such  as  the  machine  shop,  measurement  equipment, 
and  computer  facilities. 

One  immediate  necessity  was  the  fabrication  of  a galvanized  sheet  metal  and 
plywon.i  oil  trap  and  floor  protector.  Normal  operation  of  the  machine  is  messy 
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and  produces  a continual  flow  of  oil  to  4he  floor  and  surrounding  work  surfaces. 
Occassional  malfunctions  can  produce  a problematic  amount  of  oil  in  a matter  of 
the  seconds  it  takes  to  restorp  control  or  shut  off  the  system.  Connection  of  the 
hydraulic  lines  was  straight-f orward  and  quickly  accomplished.  High  pressure  hoses 
were  obtained  from  a local  supplier  by  special  order. 

Fortunately,  we  had  been  forewarned  about  the  impaling  capabilities  of  the 


machine  under  certain  extrusion  conditions  specifically,  pinch-off  or  over- 
extrusion can  results  in  a sudden  loss  of  pressure  in  the  extrusion  chamber  as 
the  end  of  the  billet  is  extruded  through  the  die.  The  next  task  was  therefore  to 
build  a suitable  safety  barrier  which  could  be  placed  in  position  to  absorb  or 
deflect  any  projectiles.  To  handle  the  design  and  construction  of  the  safety 
barrier,  we  called  upon  the  services  of  a graduate  student  on  the  project.  The 
result  of  his  work  is  a sturdy  wood  barrier  five  feet  by  eight  feet  and  four  inches 
thick  mounted  on  a two  and  a half  foot  tall  pedestal  with  casters  on  the  bottom. 

Both  the  forces  involved  in  the  operation  of  the  machine  and  the  adequacy  of  the 
barrier  were  dramatically  and  unexpectedly  illustrated  when  during  one  of  the  first 
extrusions  a loud  bang  was  followed  immediately  by  a puff  and  trail  of  smoke.  As 
the  smoke  cleared  and  the  panic  subsided,  we  observed  a two  foot  long  aluminum  rod 
quivering  dead  center  in  the  strategically  located  safety  barrier.  Figure  1.  shows  a 
photograph  of  the  hydrostatic  extrusion  machine  in  its  present  location.  The 
movable  barrier  (in  the  background  of  the  photograph)  is  in  the  proper  position 
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for  the  direct  mode  of  hydrostatic  extrusion. 

The  first  few  trial  extrusions  were  successful.  Figure  2 shows  three  sample 
extrusions  of  1100  aluminum.  In  each  case  the  original  rod  diameter  was  .750". 
Listing  from  top  to  bottom,  the  extrusion  ratios  are  16:1,  25:1  and  36:1.  Unfortu- 
nately, problems  quickly  developed  in  the  hydraulic  extrusion  piston  and  the  pressure 
intensifier.  Since  both  of  thbse  failures  were  unexpected,  it  seamed  appropriate 
to  dismantle  and  clean  the  entire  machine,  check  alignment  and  tolerances,  and 
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Figure  1.  The  above  photograph  shows  the  hydrostatic  extruder  in  its 
present  location.  Both  the  hydrostatic  control  panel  and  protection 
barrier  are  visible  in  the  background.  The  barrier  is  properly 
positioned  for  the  direct  extrusion  mode. 


i 2.  Sample  aluminum  extrusions.  Initial  billet  diameter  Was 
Listing  from  top  to  bottom,  the  extrusion  rations  are  16:1, 
and  36:1. 
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replace  all  seals.  The  original  problem  we  faced  with  the  hydraulic  piston  was 
that  the  die  stem  assembly  would  not  retract.  When  the  hydraulic  pressure  was 
boosted  to  increase  the  force  of  retraction,  the  die  stem  slipped  away  from  the. 
holding  collet  and  oil  gushed  profusely  from  both  ends  of  the  machine.  The  problem 
was  diagnosed  as  jammed  die  stem  guides  — probably  the  result  of  bent  die  stem  pins. 
Special  tools  were  machined  and  fabricated  to  facilitate  the  removal  of  the  cylinder 
internal  parts.  Disassembly  did,  in  fact,  reveal  bent  ox  scored  guide  pins.  We 
surmised  that  either  the  guides  had  not  been  aligned  properly  originally  and  that 
this  had  caused  the  pins  to  bend  as  the  piston  was  fully  retracted  or,  alternatively, 
that  the  piston  had  simply  been  over  contracted.  The  defective  set  of  pins  was  re- 
placed with  a new  set  which  was  machined  at  Stanford. 

The  pressure  intensifler was  not  recharging  properly.  It  was  not  retracting 
even  under  high  load  Both  the  intensifier  piston  and  the  cylinder  were  removed 
and  measured.  It  was  obvious  from  simply  looking  at  the  piston  and  high  pressure 
seals  that  after  only  approximately  200  cycles,  the  intensifier  was  showing  signs 
of  uneven  wear.  The  dimensions  were  definitely  already  beyond  the  tolerances 
specified  in  the  blueprints.  In  an  attempt  to  minimize  further  delays,  it  was 
decided,  in  consultation  with  Mr.  Keathley  at  Revere  Copper  and  Brass,  that  we 
would  try  to  repair  the  intensifier  at  Stanford  by  first  carefully  polishing  and 
replacing  the  worn  seals.  Both  problems  were  satisfactorily  resolved  and  the 
machine  was  returned  to  working  condition.  To  our  dismay,  the  first  attempts  to 
extrude  after  repairs  were  made  resulted  in  failure.  Typically  the  machine  would 
either  stall  at  maximum  extrusion  pressure  (250,000  psi)  or  the  material  would 
characteristically  pinch  off.  To  complicate  matters  further  the  flow  control 
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valve  was  not  operating  properly  in  the  factory  installed  position.  The  idea 
had  been  to  place  the  valve  in  parallel  with  the  extrusion  piston  so  that  the 
meter  acted  as  a bleeder  to  vent  excess  oil  to  the  tank  (the  slower  the  extrusion 


rate,  the  more  venting  to  tank).  However,  it  became  obvious  that  the  vencing 
rate  was  very  pressure  sensitive  and  rate  control  was  impossible.  The  valve 
was  originally  placed  in  parallel  to  avoid  the  loss  in  maximum  attainable  ex- 
trusion pressure  resulting  from  the  pressure  drop  through  the  flow  meter  when  it 
is  connected  in  line.  We  had  no  choice  but  to  sacrifice  a possible  five  percent 
drop  in  extrusion  pressure  for  rate  control.  The  lines  were  therefore  rerouted 
soithat  the  valve  is  now  in  series  with  the  extrusion  piston.  Rate,  although  not 
yet  entirely  steady,  does  appear  to  be  reasonably  controllable  and  pressure  inde- 
pendent. This  will  become  clearer  as  we  utilize  our  monitoring  systems. 

We  are  uncertain  as  to  why  many  extrusion  attempts  were  unsuccessful.  Cracked 
dies,  too  much  friction  between  die  and  liner  and  too  high  a reduction  ratio  for  the 
specific  extrusion  material  were  all  suspected.  After  several  failures,  close  exam- 
ination of  the  dies  supplied  by  Revere  Copper  and  Brass  revealed  longitudinal  cracks. 
Whether  these  cracks  created  the  stall  outs  or  were  caused  by  the  high  static  pres- 
sure is  not  certain.  However,  cracked  dies  can  open  up  under  pressure  and  lead  to 
a loss  of  lubricant  (that  is,  loss  of  hydrodynamic  lubrication),  the  consequence 
being  a sharp  increase  in  the  applied  load  with  an  actual  loss  in  hydrostatic  pi'es- 
sure.  Increased  friction  between  liner  and  die  would  have  the  same  effect;  that  is, 
the  applied  load  would  need  to  increase  in  proportion  to  a frictional  .increase  in 
order  for  extrusion  to  take  place.  If  the  applied  load  necessary  to  overcome  the 
wall  friction  is  beyond  the  machine's  capability,  the  machine  stalls.  For  this 
reason  die  diameter  and  work  chamber  liner  dimensions  are  checked  frequently. 

Normally  when  under  pressure  the  die  should  be  oversized  by  .0005"  to  .0010".  This 
provides  the  metal  to  metal  contact  for  hydrostatic  medium  retention,  but  still 
permits  most  of  the  applied  force  to  be  transmitted  to  the  extrusion  sample  via  the 
hydrostatic  fluid.  Friction  buildup  was  definitely  a problem  for  us  since  the  liner 
became  badly  scored.  This  necessitated  the  removal  of  the  liner  and  rehoning  to  a 
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new  inside  dimension.  Larger  dies  would  have  Co  be  machined  to  fit  or  old  ones 
copper  plated  — a practice  which  should  yield  several  extrusions  per  coating. 

In  any  event,  it  was  obvious  that  we  would  have  an  almost  immediate  need  to 
design  and  order  additional  dies  which  would  enable  us  to  extrude  to  a greater- 
variety  of  diameters.  We  decided  to  use  the  Stanford  machine  shop  facilities  and 
avoid  going  to  outside  contractors  for  this  work.  Thus  we  could  maintain  greater 
control  and  at  the  same  time  retain  flexibility  for  modifications  found  desirable. 
Once  initial  decisions  as  to  materials  and  optimum  heat  treatments  were  made,  die 
manufacture  became  routine.  In  recent  weeks  our  time  has  been  occupied  with  problems 
inherent  in  monitoring  the  several  extrusion  variables.  It  was  noted  in  an  earlier 
report  that  an  automatic  system  to  monitor  die  jacket  pressure  was  essential  for 
the  prevention  of  internal  damage  to  the  extrusion  chamber.  The  ultra  high  pressure 
developed  in  the  extrusion  chamber  can  only  be  contained  by  a corresponding  mainten- 
ance of  high  pressure  in  the  die  jacket.  If  die  jacket  pressure  were  to  drop  the 
load  chamber  liner  would  likely  split.  We  conceived  a plan  which  calls  for  mounting 
strain  gages  externally  on  the  high  pressure  cylinder  and  would  thus  correlate 
cylinder  strain  to  die  jacket  pressure.  When  a specific  lower  or  upper  limit  is 
exceeded  the  machine  should  shut  down  automatically.  This  safety  system  has  been 
designed,  but  not  built. 

For  our  initial  studies  the  most  important  problem  facing  us  is  the  accurate 
determination  of  the  extrusion  temperature.  We  would  like  to  know  if  extrusion 
conditions  are  isothermal;  and  if  not,  how  much  heating  takes  place.  Much  thought 
and  energy  went  into  devising  various  schemes  for  measuring  temperature.  The  object 
was  to  get  an  accurate  reading  at  the  die  exit  where  the  temperature  should  be  the 
greatest.  Access  to  the  die  was  a problem  since  space  is  extremely  limited.  We 
were  hoping  for  rather  high  response  times  in  order  to  monitor  possible  temperature 
fluctuations  during  stick-slip  or  strain  rate  changes  (that  is  increase  in  extrusion 
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rate  during  test).  Methods  such  as  infra-red  sensors  do  not  appear  feasible  at  the 
moment  due  to  inaccessibility  of  the  die.  The  tentative  decision  was  to  design  and 
build  a prototype  model  using  thin,  unshielded  thermocouples  (see  Figure  3).  Th< 
thermocouple  was  surrounded  by  a thermal  insulating  ring  of  teflon  to  prohibit  any 
thermal  contact  between  sensor  and  machine.  Clearance  between  sensor  and  extruded 
rod  was  to  be  no  more  than  .005".  This  close  tolerance  was  chosen  so  that  the  sensor 
would  be  in  contact  with  either  the  extruded  rod  or  the  thin  film  of  hydrostatic 
medium  which  is  also  extruded  as  a coating.  Since  monitoring  is  considerably  easier 
in  the  direct  mode  of  extrusion  where  the  die  is  stationary  our  initial  measuring 
attempts  concentrated  on  this  mode.  The  appropriate  plugs  and  dies  were  machined 
and  heat  treated  and  the  machine  set  up  to  operate  in  the  direct  extrusion  mode. 
Preliminary  tests  demonstrate  that  the  thermocouple  setup  does  work.  However,  there 
is  still  some  fluctuation  in  the  temperature.  This  is  probably  not  the  true  response 
of  the  extruded  material.  More  likely,  it  is  the  result  of  the  thermocouple  making 
and  breaking  contact  with  the  film  or  extruded  rod.  In  addition,  the  response  time 
must  be  calibrated;  and  most  importantly  the  accuracy  of  our  temperature  measurements 

must  be  determined. 

Although  the  extruder  will  be  used  mainly  as  a processing  tool,  we  want  to 
develop  the  capabilities  to  record  the  history  of  each  extrusion.  Of  course, 
knowledge  of  the  temperature  of  extrusion,  extrusion  pressure,  and  rate  will  be 
important  in  determining  the  mechanical  properties  of  these  materials.  But  this 
information  will  also  be  useful  in  analyzing  the  extrusion  process,  die  efficiency, 

lubricant  efficiency,  or  other  extrusion  phenomena. 

We  have  purchased  a linear  potentiometer  to  be  used  as  a displacement  gage  to 
position  accurately  the  extrusion  piston  versus  time  during  an  extrusion.  In  the 
same  manner  the  hydrostatic  pressure  and  extrusion  temperature  will  be  monitored 
with  respect  to  time.  Thus  a record  of  the  metal  working  history  of  each  pom'-  on 
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Figure  3.  Thermocouple  design  for  measuring  the  temperature  of  extruded 
products  at  the  die  exit. 
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the  extruded  rod  will  be  preserved  for  later  refereuee.  A Hewlett-Packard  computer 
system  with  analog  to  digital  conversion,  data  processing  capabilities,  plotters 
and  visual  display  has  recently  been  installed  in  the  Structures  Laboratory  and 
has  been  made  available  for  our  use.  He  will  be  plugging  our  strain  gages  and 
sensors  into  it  directly  for  data  storage  and  manipulation. 

Our  immediate  future  will  be  concerned  with  integration  of  the  Hewlett-Packard 
computer  into  our  extrusion  system.  Three  weeks  to  one  month  will  probably  be  re- 
quired to  set  up  the  equipment  and  learn  the  techniques  for  computer  operation  and 

commands. 

Simultaneous  with  efforts  to  implement  the  computer  analysis  will  be  trial 
extrusions  of  harder  materials.  Our  goal  is  to  develop  capabilities  for  wa,  m and 
hot  extrusions.  Battelle  reportedly  has  the  ability  to  extrude  at  1000"C  with 
pressures  to  250,000  psi.  Unfortunately,  details  of  Battelle  s technique  and  equip 
m„nt  are  proprietary.  Our  hydrostatic  extrusion  machine  was  not  designed  specifically  | 
with  hot  extrusion  in  mind  but  wo  hope  that  ideas  we  have  developed  will  permit  ex- 
trusion approaching  1000°C  and  at  pressures  of  250,000  psi.  Barring  major  unforseen 
difficulties,  extrusion  temperatures  of  500-600‘C  at  full  pressure  (250,000  psi) 
should  be  readily  achievable.  The  first  precaution  must  be  to  insure  against  over 
heating  of  the  hydraulic  oil  and  loss  of  tempered  strength  of  the  die  Jackets.  At 
higher  temperatures,  heat  retention,  thermal  shook,  thermal  expansion  and  selection 


...  1 


of  material,  hydrostatic  medium  and  lubricant  are  all  factors  which  may  be  problematic. j 
He  considered  both  internal  and  external  heating  of  the  billet.  The  alternative  . 


to  heat  externally  and  quickly  transfer  billet  and  hydrostatic  medium  offered  the 
best  chance  for  immediate  success.  Therefore,  our  efforts  have  concentrated  on 
perfecting  this  system.  Figure  4 illustrates  the  proposed  modification  of  the  pro 
sent  system  to  accomodate  warm  extrusion  (200-750"C).  He  are  presently  still  de- 


signing, checking  into  the  materials  selection  and  also  researching  other  related  j 

work  which  has  recently  been  published. 
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II . The  Basis  for  an  Elastic-Plastic  Code 
(Professor  Erastus  H.  Lee) 

Introduction 

lor  satisfactory  analytical  assessment  of  metal  forming  problems,  it 
is  necessary  to  evaluate  the  varying  stress  distribution  in  an  element 
since  forming  defects,  such  as  the  development  of  internal  cracks,  can 
depend  on  stress  history,  and  residual  stresses  can  be  important  in  deci- 
ding  the  utilization  of  a formed  part.  Elastic  characteristics  play  an 
important  role  in  the  determination  of  stress,  even  in  combination  with 
extensive  plastic  flow  which  may  involve  strains  a thousand  times  elas- 
tic strain  magnitudes.  Thus  analysis  of  metal  forming  problems  for 
such  assessments  must  be  based  on  elastic-plastic  theory.  The  same 
is  true  for  other  stress  analysis  problems  when  plastic  flow 
occurs. 

Because  the  plasticity  laws  are  incremental  in  nature,  they  result 
in  relations  between  stress-rate  and  strain-rate,  or  equivalently  in 
numerical  evaluations,  between  stress  and  strain  increments.  For  the 
commonly  adequate  rate-independent  laws,  linear  relations  between  stress- 
rate  and  strain-rate  arise,  the  coefficients  being  functions  of  the  cur- 
rent stress  for  the  common  laws  when  plastic  flow  is  taking  place,  and 
otherwise  the  elastic  laws  apply  in  incremental  form.  Because  of  the 
structure  of  these  laws,  elastic-plastic  problems  are  commonly  solved  in 
terms  of  equations  for  stress-rates  and  strain-rates,  containing  stresses 
as  coefficients.  Thus  a time  step  forward.  At  , from  the  current  situa- 
tion at  time  t gives  the  solution  at  time  t + At  with  stress  o(t)  •!- 
oAt  , and  similarly  for  other  variables.  Then  a new  time  step  can  be 

o 

taken  and  the  process  repeated,  a is  the  appropriate  stress-rate. 
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Extensive  studies  of  the  application  of  these  laws  to  stability  and  unique- 

•k 

ness  of  solutions  have  been  made  by  Hill  (see,  for  example  [1],  [2],  [3]  ) 
where  he  shows  that  care  in  the  selection  of  stress  definitions  and  stress- 

rate  and  strain-rale  expressions  is  important  for  a satisfactory  develop- 

> 

raent  of  the  theory.  Rice  [4]  has  pointed  out  that  such  questions  are  also 
important  in  developing  a satisfactory  theoretical  basis  for  elastic-plastic, 
stress  analysis,  particularly  in  the  common  circumstance  that  the  tangent 
modulus  in  plastic  flow  is  of  the  order  of  the  stress.  Conveeted  and  rota- 
tion terns  then  become  important  in  the  stress-rate  expression,  and  analo- 
gously stress  variables  should  be  selected  so  that  the  influence  of  rate 
of  deformation  of  the  boundaries  of  the  body  does  not  affect  the  variational 
principle  which  replaces  the  equilibrium  equation.  This  requirement  can  be 
achieved  by  using  the  unsymractric  nominal  stress  (Piola-Kirchoff  I)  in  which 
the  stress  is  defined  as  force  per  unit  undeformed  area.  The  variational 
principle  then  involves  an  integral  over  the  undeforraed  body  which  is  fixed. 

Plastic  flow  is  essentially  a fluid  type  phenomenon  which  can  be  most, 
conveniently  expressed  in  terms  of  the  current  conf iguti._ion  of  the  material. 
Thus  a reference  configuration  which  remains  invariant  throug'ntout.  the 
motion  is  not  appropriate  and  so  the  current  configuration  is  adopted  as 
the  reference  state  for  evaluation  of  the  deformation  from  t to  t + At  , 
where  At  is  sufficiently  small  for  adequacy  of  first  order  theory. 

The  framework  described  above  provides  a satisfactory  foundation  for 
a finite-element  clastic-plastic  code  as  discussed  by  McMeeking  and  Rice  [5). 
In  effect,  by  choosing  the  current  configuration  as  the  reference  state  the 


Numbers  in  square  brackets  refer  to  the  bibliography. 


Cauchy  stre.ua  (or  true  stress  in  Cartesian  coordinates)  the  u, asymmetric 
nominal  stress  (Lagrange  or  Piola-Kirchhoff  1)  and  the  symmetric  nominal 
stress  (KirchhoH  or  Piola-Kirchotf  II)  all  have  identical  values  at  the 
current  time  mhich  facilitates  utilisation  of  the  appropriate  stress  for  | 

the  appropriate  component  of  the  calculation.  Although  the  stress  com-  » 
ponents  themselves  are  identical,  rates  of  change  of  the  different  stresses 

are  not  the  same. 

Development  of  the  Theory 

Following  Hill  [2]  and  using  for  the  most  part  his  notation,  we  con- 
sider the  unsymmetric  nominal  stress  (variously  referred  to  as  Lagrange  or 
Piola-Kir chof £ I)  s defined  so  that  the  jth  component  of  the  iorce 
transmitted  across  a deformed  element,  which  in  the  initial  or  reference  | 


O 

state  had  area  dS  and  unit  normal  > 1S 


dS  v.  s . ~ dF 
i ij  J 


Hill  considers  (P.  214  of  [2])  rate  or  flow  type  constitutive  laws  of  the 


' _ __9E (2) 

’ij  30V./3X.) 
j .1  i 

vrhere  X are  rectangular  Cartesian  coordinates  in  the  initial  or  refer- 
ence configuration,  E is  a homogeneous  function  of  degree  two  in  the 
velocity  gradients,  8v./3X.  , and  where,  s±.  is  the  partial  time 
derivative  at  fixed  X , i.e.  a material  derivative  at  a particle.  The 
velocity  v^(X,t)  gives  the  distribution  at  time  t expressed  xn  the 
initial  coordinates  of  the  corresponding  material,  points  (note  that  a 
tilde  under  a symbol  denotes  vector  or  tensor  in  absolute  notation). 


Wc  consider  boundary  value  problems  in  which  for  a volume  V in  the 


reference  state,  at  time  t stress  rates 


sought  for  prescribed  nominal  traction  rates  F.  over  the  part  of  the 


velocity  v.  over  the  remainder  of  the  surface  S 


anu 


g.  per  unit  initial  volume.  Then  the  variational 


principle 


in  the  class  of  continuous  differentiable  velocity  fields  satisfying 

O 

the  velocity  boundary  condition  on  S , characterizes  the  solution, 
for  it  yields  the  equilibrium  equations 


and  the  boundary  traction  rate  condition 


for  nominal  stress,  s_„  , and  the  reference  geometry. 

In  writing  the  elastic-plastic  constitutive  relation,  we  wish  to 
associate  the  velocity  gradient  in  (2)  with  the  rate  of  deformation  or 
velocity  strain: 


where  x.  are  Cartesian  coordinate; 


expressing  the  position  of 


'31 
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particles  in  the  deformed  body 


Xi  13  Xi^X,t^ 


(7) 


Thus  to  permit  simultaneous  use  of  (3)  and  the  plasticity  laws  expressed 
in  the  usual  form  of  rate  of  defamation  of  the  current  configuration,  •#. 
Hill  takes  the  current  configuration  to  be  the  reference  state,  and  hence 


\ = x.(X,t) 


(8) 


for  a particular  t . The  theory  is  expressed  in  this  form  for  evalua-- 

tion  of  s..  and  v.  and  hence  the  solution  and  configuration  at  t + At 
13  1 

which  provides  a new  reference  state  for  evaluation  of  the  next  time 
step  At  . 

Note  that  at  the  instant  t , when  the  current  and  reference  con- 
figurations are  identical,  the  nominal  stress  components  s _ are  equal 
to  the  Cauchy  or  true  stress  components  » so  that  at  this  instant 


’ij 


is  symmetric. 


The  device  of  selecting  the  configuration  at  time  t to  be  the 
reference  state  for  evaluation  of  the  solution  at  time  t + At  thus  per- 
mits simultaneous  use  of  the  convenient  variational  principle  (3)  in 

terms  of  nominal  stress  s..  and  a fixed  geometry  and  the  familiar 

ij 

plasticity  laws  expressed  in  teirms  of  the  Cauchy,  or  true  stress,  a . . 

By  working  with  curvilinear  convected  coordinates,  , having  an 

arbitrary  configuration  in  the  reference  or  initial  state,  Hill ([2], 
p.  219  ff.)  shows  that  the  rate  potential  (2)  generates  associated  rate 
potentials  for  other  stresses  and  stress  rate  expressions,  some  of  which 
are  more  convenient  for  representing  elastic  plastic  laws. 


Consider  curvilinear  coordinates  E,  in  the  initial  or  reference 

o 

state  with  base  vectors  e . After  deformation  as  shown  in  Fig.  1 

°a 

a' 

these  become  £ , having  the  same  values  for  the  same  material  particles 

* 

as  . Then  £G  are  the  convected  coordinates  and  the  corresponding 


base  vectors  are  g , , which  are  g^  deformed  by  the  motion.  Then  the  *. 


a 


Cauchy  stress  tensor  a has  contravariant  components  a ” in  convected 
coordirates  such  that  the  force  dF  transmitted  across  an  element  of  the 
deformed  body  of  area  dS  and  unit  normal  vector  v * is  given  by: 


.dF  = oa  ^ (va-'dS')gg^ 


(9) 


Since  the  primed  coordinates  are  made  evident  by  the  notation  v and  g 
for  normal  and  base  vectors  in  the  deformed  state,  the  primes  will  usually 


at  ct  . 

be  dropped  hereafter,  for  E,  = £ denote  the  same  material  point. 


For  the  Lagrange  or  Piola-Kirchhoff  I stress,  s^.  , the  force  com- 


puted in  the  reference  frame  is  that  actually  acting  across  the  deformed 
element  (see  Fung  [6]  p.  437  for  the  usual  definition  in  Cartesian  coordi- 
nates). For  the  present  consideration  of  convected  coordinates,  the 
expression  for  dF  thus  takes  the  form 


aS  0 00 

dF  = s“P(v  dS) g 

•*>  Ct  *w  P 


do) 


For  the  Kirchhoff  or  Piola-Kirchhoff  II  stress,  x . . , Fung  points  out. 

ij 


that  the  force  vector  computed  in  the  reference  frame  must  be  transformed 
by  the  motion  to  give  the  actual  force  across  the  deformed  element,  so 
that  for  convected  coordinates 


O n O O O 

dF  = x P(v  dS)Bft 
v ct  7 “(5 


17, 
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becomes 


dF  = xa^(vudS)gp  (11) 

Nanson's  relation  for  area  element: 

c<  o o 

pv  dS  = pv  dS  (12) 

a 1 a 

then  gives,  using  (9)  and  (11) 

= l aaf!  - J o'"3  (13) 

P 

where  J is  the  Jacobian  of  the  transformation  from  the  reference  to 
the  deformed  state,  i.e. 


X ?y  “ J^a 


O 

* 8 J 


(14) 


Fig.  2 shows  the  particular  situation  when  the  5°  coordinates  in 
the  reference  frame  are  Cartesian,  Xa.  The  Cartesian  coordinates  repre- 
senting points  in  the  deformed  body  according  to  the  point  transformation, 
equ.  (7),  are  x . The  usual  definition  of  the  Kirchhoff  or  Piola-Kirchhoff  II 
stress  (see  Fung  [6]  p.  439)  is  given  in  terms  of  this  point  transformation. 
From  equ.  (7),  define  the  deformation  gradient 


F 


(15) 


Then  the  Piola-Kirchhoff  II  stress  t is  given  in  terms  of  the  Cauchy 
stress  in  Cartesian  coordinates  x , a by 


t » J f"1  o f"1  ([7]  p.  125)  (16) 

M <N>  »J 

a8 

where  J = ^iet  (F) „ This  is  in  accord  with  (13)  where'  o are  the 
contravar lant  components  of  the  Cauchy  stress  with  respect  to  the 


convected  coordinates  X , for  if  a are  the  Cartesian  coordinates 
of  o with  respect  to  x , then  the  tensor  change  of  variables  law  gives 


3x“'  axB'  SiJ 

ox-  3x3 


in  terms  of  the  coordinate  transformation  from  x X'  in  the  deformed 

a a' 

geometry.  Because  of  the  property  of  connected  coordinates  that  X = X 
for  the  same  particle,  (7)  also  expresses  the  coordinate  transformation, 
and  (13)  and  (17)  are  seen  to  be  equivalent  to  (16) . Note  that  s ' and 
x°^  are  tensor  densities  and  not  absolute  tensors,  so  that  equations  such 
as  (2)  are  not  pure  tensor  relations 

For  rate  independent  constitutive  relations  the  rate  potential  func- 
tion E in  (2)  is  a homogeneous  function  of  degree  two.-  For  the  choice 
that  the  reference  state  is  the  current  state,  (2)  follows  from  the  exis- 
tence of  an  associated  rate  potential  function  F(cQp)  , of  the  strain-rate 
components 


r — (v  + v ) 

a 3 2 v «,3  P,ay 


where  v is  the  velocity  vector  in  the  deforming  body  and  the  comma  denotes 

• aB 

govariant  differentiation.  This  function  generates  x through 


*a3  3F 
T = 3e 


where  the  superposed  dot  indicates  time  derivative  of  the  convected  com- 
ponents or  convected  derivative.  This  is  equivalent  to  the  partial  time 
derivative  at  fixed  X , or  material  derivative.  The  structure  of  (19) 


indicates  that  since  x is  a tensor  density,  F is  a scalar  density  and 


19. 


rot  an  absolute  scalar  invariant.  Some  details  of  the  derivations  needed 


are  given  in  the  Appendix,  in  particular  for  the  relationship  between 
‘aS  . *afJ  , . , , , . 

s and  T which  for  the  particular  choice  of  reference  state  mentioned 


ora 


multiplying  bouh  sides  by  v 


[3F/3(e  ) ]e  0 + tutv 


using  (18)  and  the  fact  that  t is  symmetric.  Now  F(e  „)  is  also 

ap 

a homogeneous  function  of  degree  two  since  for  plasticity  (19)  is  rate 
Independent,  hence  Euler's  theorem  for  homogeneous  functions  permits 
us  to  write  (21)  in  the  form 


Again  using  Euler's  theorem,  this  is  consistent  with  (2)  in  terms  of 


convected  coordinates 


With  t 


k town,  this  establishes  that 


a homogeneous  function 


econd  degree  in  v 


The  laws  of  plasticity  arc  normally  obtained  by  measuring  true 
stte88»  and  increments  of  strain  defined  in  terms  of  Cartesian  coordinates 
in  the  current  conf iguration  without  rotations  occurring.  Since  the  theory 
must  apply  in  the  presence  of  rotations,  their  influence  must  not  affect 

the  stress  rate  term  in  the  constitutive  relation,  hence  a Jaumann  or 

oAnmnriTfe  Uc  will  work  in  tei'ms  of 
spin-invariant  .'.tress  rate  seems  appropriate,  we  win 

. . . t.  a-i  rh  -Is  the  reference  eonf iguration,  and 

the  configuration  at  time  t , w.n  ch  is  tne  reiei 

utilize  Cartesian  coordinates  x . ^ . 

in  formulating  the  finite  element  theory  for  numerical  analysis  of 

elastic-plastic  problems,  we  wish  to  use  the  variational  relation  (3)  xn 
terms  of  s..  and  a constitutive  relation  associated  with  C 1=0  terms 
of  T since  it  can  be  conveniently  associated  with  measurements  of  plas- 
ticity'lows.  Wo  have  seen  that  a law  in  the  form  (19)  implies  the  validity 
of  (2)  and  hence  the  variational  principle.  As  pointed  out  in  15)  this 
structure  in  terms  of  T.  . leads  to  symmetric  stiffness  matrices  in  the 
finite  element  formulation  whieh  simplifies  the  numerical  procedures. 

Now  the  relationship  between  the  nominal  stress  s and  the  Cartesian 
e . 

true  stress  a IB 


F s « J o 

~ •»  ' 

see,  for  example  [7],  p.  125,  where  the  nominal  stress  is  definedas  the 
transpose  of  s (or  it  eon  be  deduced  from  (9),  (10)  and  (17)).  Taking 
the  material  derivative  of  (24),  and  noting  that  1 and  J ore  unity 
coineldent  reference  and  current  configurations,  one  obtains 


The  difference  between  the  Jaumann  derivative  of  a and  its  material 
derivative  is  the  contribution  of  the  rotation  of  the  axes  which  rotate 


with  the  body  according  to  the  anti-synunet:ric  tensor 


to  denote  the 


where  A denotes  the  anti-symmetric  part.  Using 


Jaumann  derivative,  this  determines 


Combination  of  (25) 


Strictly  speaking,  t only  defines  the  Kirchhoff  stress  for  convectea 
coordinates.  Being  a tensor  density,  definition  for  other  coordinates 
does  not  follow  the  tensor  relation  (17)  which  yields  other  components 
of  the  stress  tensor  a . However,  we  wish  to  define  the  first  term  on 


/Dt  where 


and  this  can  be  written 


This  result  follows  by  talcing  the  Jaumann  derivative 


the  product  and  noting  that  the  derivative  of  the  scalar  density  J is 
unaffected  by  the  rotation  of  the  axes  and  that  the  instantaneous  value  oi 
J is  unity.  Note  that  to  define  a Jaumann  derivative  of  x , J must  be 
defined  for  other  than  convected  coordinates  since  rotated  rigid  coordinates 


can  be  written 


are  involved 


How  combining  (20)  and  (29) 


In  view  of  (19)  written  for  initial  Cartesian  coordinate 


with 


homogeneous  of  second  order  in  the  strain  rates,  multiplication 


of  (30)  by  D. . and  using  Euler 


theorem  gives 


Thus  defining 


and  using  Euler's  theorem,  a third  rate  potential  function  is  generate' 


The  classical  elastic-plastic  isotropically  work  hardening  law 


(Prandtl-Reuss) , commonly  giving  strain-rate  as  a linear  function  of 


stress-rate,  can  be.  inverted  to  give  stress-rate  and  takes  the  form 


where  E and  v are  Young's  modulus  and  Poisson's  Ratio,  o'  denotes 


-v-  ii1 


«y.  -4,  ‘^' 


stress  deviator,  o is  the  current  tensile  yield  stress  and  h the  current, 
gradient  of  the  true-stress  logarithmic  plastic  strain  curve  in  a tension 
test.  The  Jaumann  derivative  is  used  for  stress-rate  as  mentioned  above  in 
order  to  eliminate  rotation  effects,  and  the  last  term  in  the  brackets  is 
dropped  when  the  increment  of  deformation  is  clastic.  We  have  shown  that*. 


there  is  a rate-potential  function  G for  i _ , 


but  one  does  not  exist 


for  a . This  means  [5]  that  a non-symraetric  finite-clement  stiffness 

ij 


matrix  would  be  deduced  using  (34).  Replacing  o_^_,  by  ia  the  stre.s^ 


rate  term  would  yield  a symmetric  stiffness  matrix  which  would  simplify  the 
numerical  analysis.  Moreover  such  a change  is  appropriate  in  terms  of  its 
representation  of  the  physical  laws.  The  J term  in 


ij 


= J a 


(35) 


.ij 


arises  in  more  accurate  representation  of  the  laws  of  elasticity  than 
Hooke's  law  used  in  classical  elasticity.  It  is  associated  with  geometri- 
cal non-linearity  which  expresses  the  non-linear  influence  of  finite  strain. 
Xn  effect  it  expresses  the  fact  that  energy  density  should  be  referred  to 
a unit  initial  volume,  for  a unit  current  volume  implies  change  of  energy 
density  simply  because  the  amount  of  material  containing  -it  changes.  For 
example,  the  term  (p/pQ)  appearing  in  equation  (26)  of  [9]  is  equivalent 
to  replacing  a by  Ja  , and  it  wa s pointed  out  in  that  paper  (p.  935) 
that  such  a term  provides  a good  approximation  to  non-linear  elasticity 
metals  with  input  of  only  the  two  classical  elastic  constants.  A similar 
modification  of  the  classical  laws  of  plasticity  was  suggested  in  [10] 
where  it  was  found  compelling  to  express  the  yield  stress  in  terms  of  Jo 
(equs.  (33)  and  (34)  of  [10]).  Again  this  was  based  on  the  requirements 


1 


m 


* 


i 


. - .....  » • - * j 

- ..  — x,  ^ — — - - 

■ ..... ...  ..  •• ...  • ...  ...  . .... 


of  geometrical  non-linearity,  which  of  course  are  independent  of  specific 
material  characteristics. 

Introduction  of  (2)  into  the  variational  principle  (3)  expresses 
it  in  the  form 


3v.  o , o f „ 0 

/ s . Sfc1)  d V - / F v dS  - / g v dV  = 0 

i i|  o .1  .1  t>  j j 


vihich,  since  the  reference  state  is  the  current  state  (X  = x)  can  be 


v?ri  tten 


3v„ 

/ s 5 (~i)dV 
V « Jxi 


/ F.v.dS  - / g.v.dV  - 0 


Substitution  of  (29)  and  applying  algebraic  manipulation  based  on  sym- 
metries then  yields  the  variational  principle  in  the  form  (see  equ.  (5) 
of  [5]) 


/ 5<V 


- - o . .6  (2D#1  D ” v v,  . ]dV 

2 ij  lit  kj  k,i  k,j 

~ / F.v.dS  - / g.v.dV  = 0 

* -I  l * 1 T 


where  the  subscript  ,i  denotes  the  operation  * ^ow  utilizing 

(33),  the  principle  takes  the  form  (see  equ.  (15)  of  [5]). 
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By  (33),  and  using  Euler's  theorem 
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(AO) 


since  3G/3D  is  homogeneous  of  first  degree  in  D , hence  £..]cl  is 
symmetric  in  ij  kl  as  well  as  i j and  k 1 . Since  G is 


homogeneous  of  second  degree 


1 ®?ij 

G(D)  - f s?  D 


(41) 


The  symmetry  implicit  in  the  variational  principles  (38)  and  (39)  with 
(AO)  and  (Al)  imply  symmetric  stiffness  matrices  which  carry  through  to 
the  finite-element  formulation  [5], 


Discussion 


As  discussed  by  Rice  [A]  and  McMeeking  and  Rice  [5],  the  development 
reviewed  in  these  notes  brings  out  the  importance  of  convection  effects 
and  w propriate  stress-rate  definitions  in  formulating  elastic-plastic 
theory.  Thus  the.  precise  analysis  of  continuum  theory  must  be  applied  to 
obtain  reliable  results.  Common  snail  strain  assumptions  are  not  valid 
even  for  incremental  theory  based  on  small  strain  increments.  Ihis  situ- 
ation is  implicit  in  comparison  of  the  first  and  second  terms  in  tne  first 
volume  integral  in  the  variational  principle  (38) . In  plastic  flow  the 
coefficient  of  SD  in  the  first  term  is  ~ h D , where  h is  the  gradient 
of  the  tensile  stress-plastic  strain  curve,  whereas  the  second  term  is 
~ oD6D  . The  second  term  and  the  difference  between  the  Jaumann  derivative 
and  other  simpler  time  derivatives  can  only  be  neglected  if  h » ci  . 
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The  relative  error  in  neglecting  such  terms  is  effectively  independent  oi 
the  magnitude  of  the  strain  increment  adopted,  so  that  small  steps  do  not 
permit  simplification  in  this  regard.  For  many  metals  h ~ a • 

It  is  interesting  to  note  that  the  complications  which  arise  in  elastic 
plastic  analysis  result  from  the  elastic  component  of  strain,  and  not  the" 
Plastic.  For  stress  and  strain  deviators,  elastic-ideally  plastic  defor- 


mation with  a Mises  yield  condition,  J 


where  X is  a parameter.  Multiplication  of 


gives 


For  isotropic  work  hardening  with  a Mises  yield  condition, 


The  first  terms  on  the  right  hand  sides  of  (42)  and  (44)  are  the  elastic 
strain  rate  components.  Thus  rigid-ideally  plastic  theory  ((42)  with  the 
o term  deleted)  gives  a relation  between  stress  and  velocity  gradient 
with  no  complications  due  to  stress  rate,  which  greatly  simplifies  the 
analysis,  apart  from  the  difficulty  of  determining  the  rigid  regions. 

o 

Similarly  for  work  hardening  rigid-plastic  analysis  ((44)  with  the  a term 
deleted),  only  the  rate  of  change  of  a stress  invariant  occurs,  which  is 
simpler  to  include  than  & tensor  rate.  For  elastic-plastic  theory  it  is 
the  elastic  term  which  introduces  stress-rate  and  the  consequent  complication 


Many  technologically  important  metal-forming  problems  are  steady  state 


In  planning  to  use  an  analysis  of  the 


ituations,  it  is  fortunate  that  the 


type  considered  here  to  evaluate  such 


stress-rate  terms  appearing  in  the  variational  principle 


teady  case,  hencu  singular  computational  condl 


do  not  approach  aero  in  the 


not  the  case  in  some  simpler  and 


tions  need  not  be  anticipated.  This  i 


inadequate  approaches  to  this  problem  in  which  sufficient  care  wa 


devoted  to  the  appropriate  choice  of  stress-rate  definition 
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Appendix 


For  simplicity, relations  are  developed  for  Cartesian  axe 


reference  state  which  is  instantaneously  coincident  with  the  current 
state.  The  general  theory  for  convected  coordinates  carries  through 


a similar  manner  but  can  be  technically  much  more  involved 


Equations  (16)  and  (24)  yield  the  relation 


Material  differentiation  of  this  relation,  which  is  equivalent  to  con 
vected  differentiation,  gives 


and  for  the  special  reference  configuration 


hence 
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